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Summary. We examine effects of small-scale fluctuations with angle in the neutrino 
radiation in core-collapse supernova explosions. As the mode number of fluctuations 
increases, the results approach those of spherical explosion. We conclude that global 
anisotropy of the neutrino radiation is the most effective mechanism of increasing 
the explosion energy when the total neutrino luminosity is given. 



1 Introduction 

o 

ft has been turned out that most simulations of core-collapse supernova explosions 
' with spherical symmetry fail to produce a successful explosion[l]. In addition, ob- 

f"^ , servations suggest that the ejecta of core-collapse supernova explosions are aspher- 

ical(e.g., [2]). These facts lead us to multidimensional simulations. In 2-D and 3-D 
simulations [3, 4, 5, 6, 7, 8, 9, 10, 11], it has been shown that multidimensional effects, 
such as convection inside the proto-neutron star and convective overturn around the 
O ' neutrino-heated region, increase the explosion energy and can trigger a successful 

explosion [4, 6, 12]. 

When a proto-neutron star rotates, the neutrino flux is expected to be enhanced 
along the rotational (polar) axis. Janka and M6nchmeyer[13, 14] first discussed the 
possibility of aspherical neutrino emission from a rapidly rotating inner core. They 
argued that a neutrino flux along the polar axis might become three times greater 
r> ' than that on the equatorial plane. 

Shimizu et al.[15, 9] proposed that the anisotropic neutrino radiation should 
play a crucial role in the explosion mechanism itself. They carefully investigated the 
effects of anisotropic neutrino radiation on the explosion energy. They found that 
only a few percent enhancement in the neutrino emission along the pole is sufficient 
to increase the explosion energy by a large factor, and leads to a successful explosion. 
They also found that this effect saturates around a certain degree of anisotropy. It 
should be noted here that the assumed rotational velocity of the inner core is very 
different between Janka et al.[13, 14] and Shimizu et al.[15, 9]. 

In the work of Shimizu et al.[9], they considered only a form of global anisotropy; 
the maximum peak in the neutrino flux distribution was located at the pole and the 
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minimum at the equatorial plane. On the other hand, Burrows et al.[5] have sug- 
gested that the neutrino flux can fluctuate with angle and time due to gravitational 
oscillation on the surface of the proto-neutron star. In this work, we introduce such 
small-scale fluctuations in the neutrino flux in our numerical code by modifying the 
angular distribution of the neutrino flux. We aim to study the effects of these small- 
scale fluctuations on the shock position, the explosion energy, and the asymmetric 
explosion. Details are found in ref.[16]. 

2 Numerical Code 

We perform 2-D simulations by solving hydrodynamic equations in spherical coordi- 
nates. A generalized Roe's method is employed to solve the hydrodynamic equations 
with general equations of state (EOSs). The details of our numerical technique, to- 
gether with the EOS and the initial condition used, are described in the previous 
article [9]. In our study, we have improved the numerical code of Shimizu et al.[9]; the 
cells in the ^-direction were shifted by half of the cell size[17] in order to avoid a nu- 
merical error near the pole, although the error was not serious for the investigation 
of the explosion energy. 

In the present work, the local neutrino flux is assumed to be[16] 



where a is the Boltzmann constant, and T v is the temperature on the neutrinosphere. 
In equation (1), ng represents the number of waves in the ^-direction. The case of 
no — 1 corresponds to the global anisotropy, namely, no fluctuation. C2 is a parameter 
which is related to the degree of anisotropy in the neutrino radiation. We see in 
equation (1) that the neutrino fluxes in the x (equatorial) and z (polar) directions 
become l x = l v (r,8 = 90°) oc ci and l z = l v (r,8 = 0°) oc ci(l + C2), respectively. 
The degree of anisotropy l z /l x is then represented as 



The value of ci is calculated from C2 and ne so as to adjust the total neutrino flux 
to that in the spherical model at the same T v . 

It should be noted here that the amplitude of fluctuations in the neutrino flux 
distribution for an observer far from the neutrinosphere and that on the neutrino- 
emitting surface are different. When we observe the neutrino flux far from the neu- 
trinosphere, the local neutrino flux is seen as equation (1). On the other hand, the 
neutrino flux on the neutrino-emitting surface has a profile similar to equation (1) 
but the amplitude is different. In the latter, C2 is replaced by a where a is a parameter 
which represents the degree of anisotropy of the neutrino flux on the neutrinosphere. 
It is preferable that we compare the results for the same value of a, since a is more 
directly related to explosion dynamics. The value of C2, therefore, is calculated from 
a given a, depending on ne. Although it is difficult to calculate the exact relationship 
between C2 and o, we can estimate it by assuming that the strength of the neutrino 
flux on the neutrinosphere is approximated by a profile of step function. This makes 
it possible to relate the value of C2 to the value of a for each ng. For detail, see 
Madokoro et al.[16]. 
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— = 1 + C 2 . 
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We set a — 0.31 in this work. This value of a is chosen in such a way that the 
value of l z /l x for the global model (n$ = 1) becomes 1.10. The values of C2 for each 
fluctuation model (n$ = 3, 5) are accordingly calculated. These are summarized in 
Table 1. The neutrino temperature on the neutrino-emitting surface T v is assumed 
to be 4.70 MeV. In our simulation, we have 500 nonuniform radial zones which cover 
from 50 to 10000 km in radius. For 9, we have 62 uniform angular zones from = 0° 
to 6 = 90° with equatorial symmetry for ng = 1 and 3. Note that we use a double 
value of the angular zones for n$ = 5, which is different from that used in ref.[16]. 



3 Results 

Figure 1 depicts the contour maps of the dimensionless entropy[9] distribution with 
the velocity fields for the three models at t ~ 250 ms after the shock stall. The shock 
front is represented by the crowded contour lines at r ~ 2000—3000 km for the model 
A1-T470, r ~ 1600 -2200km for the model A3-T470, and r ~ 1400km for the model 
A5-T470. We see in Figure 1 that the shock front is largely distorted in a prolate 
form for the model A1-T470. This is because the neutrino heating along the pole is 
more intensive than that on the equatorial plane. Due to increased pressure in the 
locally heated matter near the polar axis, the shock front along the pole is pushed 
up, resulting in a prolate deformation. We find that the shock position becomes 
less extended than that of the global anisotropy, and the shock front approaches 
the spherical shape when the mode number of fluctuation increases. This trend is 
especially remarkable for the model A5-T470 in which the shock front is almost 
spherical. 

Figure 2 shows the evolution of the explosion energy, as well as the thermal, 
kinetic, and gravitational energies for the three models. The energy gain for the case 
of no = 1 is the highest among others at all stages of the explosion. It is also seen 
in Figure 2 that the explosion energy decreases as the mode number of fluctuations 
in the neutrino flux increases and finally approaches that of the spherical explosion. 

Thus, we found that there are remarkable differences in the explosion energy 
depending on the mode number of the fluctuations. We also found that larger num- 
ber of modes in the fluctuations makes the result closer to that of the spherical 
explosion. This is because any small-scale fluctuations on the neutrinosphere are 
greatly averaged out when the neutrino emission is observed far enough from the 
neutrino-emitting surface. Moreover, we found that a certain broad space is needed 
to be heated by neutrinos to revive the stalled shock wave rigorously and hence the 
global anisotropy (no = 1) is the most effective to increase the explosion energy. 
These results support the claim by Shimizu et al.[9]. 



4 Conclusion 

We have investigated the effects of small-scale fluctuations in the neutrino flux on 
the core-collapse supernova explosion. The profile of the neutrino radiation field was 
specified taking its geometric effects into account. Since the small-scale fluctuations 
are averaged out for radiative and hydrodynamic reasons, the results of the fluc- 
tuation models become closer to that of the spherical explosion. Consequently, the 
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global anisotropy is the most effective mechanism in increasing the explosion energy 
when L„ is given. This supports the claim made by Shimizu et al.[9]. 



Model 


rig 


C2 


A1-T470 


1 


0.100 


A3-T470 


3 


0.051 


A5-T470 


5 


0.035 



Table 1. Simulated models. 
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Fig. 1. Contour maps of the dimensionless entropy distribution and the velocity 
fields for the three models; top left: n$ = 1 (model A1-T470) at t = 244ms, top 
right: n g = 3 (model A3-T470) at t = 254ms, bottom: n g = 5 (model A5-T470) at 
t = 252ms. 
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Fig. 2. Evolution of thermal and kinetic energy (E t h + Skin), gravitational energy 
(£ g rav) and explosion energy (E exp {) for the three models; solid line: ne = 1, short- 
dashed line: = 3, long-dashed line: ne = 5. For comparison, the result of the 
spherical explosion at the same T v is also plotted (dotted line). 
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